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Recurrent Acceleration in Axion Dilaton Cosmology 
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A class of Einstein-dilaton-axion models is found for which almost all flat expanding homogeneous 
and isotropic universes undergo recurrent periods of acceleration. We also extend recent results on 
eternally accelerating open universes. 
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I. INTRODUCTION 

The discovery that the expansion of our Universe is 
accelerating, presumably due to "dark energy" , together 
with the compatibility of recent astronomical observa- 
tions with the hypothesis of an earlier inflationary epoch, 
has given a new impetus to studies of cosmologies driven 
by scalar fields. Models with a single scalar field have 
a long history. In more recent times, models with two 
or more scalars have attracted much attention; in par- 
ticular models with both dilaton and axion [3. Of most 
relevance here is the work of Billyard et al. [2| because 
the models we study include one studied by them, and 
our methods are similar. 

We suppose that a dilaton field a and axion field Xi 
which together define a map from a £>-dimensional space- 
time to a hyperbolic space with an 57(2; R)-invariant 
metric, are coupled to gravity. The S7(2; M)-invariance 
must be broken by any non-constant scalar field poten- 
tial, but if we require only invariance under dilatations 
then an exponential potential for the dilaton is allowed. 
These considerations lead to a low-energy Lagrangian 
density of the form 



C = \J — det g 



-Xo 



where // and A are constants, as is the 'cosmological con- 
stant' A, which we here assume to be positive. We may 
choose A > without loss of generality but then the sign 
of /i is significant, so we allow for either sign. Note that 
\n\ is proportional to the radius of the hyperbolic target 
space. When /x = 0, the target space is flat and the axion 
field decouples; the resulting cosmologies were studied in 
[3, 0, [f| so we assume here that \i ^ 0. Models of this 
type can be motivated in various ways. For example the 
Freedman-Schwarz D — 4 supergravity theory [6| has a 
consistent truncation to p.l[) with 2A = — fi = 2. 

We shall study flat FLRW cosmological solutions of 
the above model as a function of the parameters (i and A. 
The D = 4 model with A = — fi = 2 was previously stud- 
ied by Billyard et al. with a string-theory motivation 
that led them to consider a 'string-frame' metric rather 
than the Einstein- frame metric implicit in (jl.ip . They 
found an interesting quasi-cyclic evolution of flat cosmo- 
logical models that they interpreted as successive peri- 
ods of expansion and contraction. As we point out here, 



the Einstein-frame universe is either eternally expanding 
or eternally contracting, but each undergoes successive 
periods of acceleration. The 'string-frame' cosmologies 
of 0] also undergo successive periods of acceleration be- 
cause this is an immediate consequence of re-expansion 
following a period of contraction, but acceleration in non- 
Einstein frames does not imply acceleration in the Ein- 
stein frame; in particular, a flat Einstein-frame universe 
cannot smoothly pass from contraction to expansion (see 
e.g. 0) so this source of cosmic acceleration is not avail- 
able in the Einstein frame. We shall show that, in Ein- 
stein frame, almost all flat expanding FLRW cosmologies 
undergo quasi-cyclic alternation between periods of ac- 
celeration and deceleration provided that both A and —ji 
are sufficiently large. The late-time behaviour is a power- 
law expansion which may be accelerating or decelerating, 
depending on the value of n/X. 

Of course, there is no physical significance to the choice 
of frame within the models we consider but this choice be- 
comes significant once one considers the additional fields 
needed for a realistic theory. In this larger context, the 
choice of Einstein frame is equivalent to the assumption 
that the energy-momentum-stress tensor of the matter 
fields responsible for structure couple to the Einstein- 
frame metric, in the sense that the motion of a free par- 
ticle associated with such a field is geodesic in this frame, 
an assumption that is consistent with observations. 

The model defined by jTT]) can be consistently trun- 
cated by taking x — Xo, for any constant xo- Solutions of 
the resulting Einstein-dilaton model are therefore also so- 
lutions of the Einstein-dilaton-axion model, and this pro- 
vides a useful check on our results. The Einstein-dilaton 
model has been much studied in the past, but some fea- 
tures in relation to eternal cosmic acceleration of open 
FLRW universes in models arising from hyperbolic com- 
pactification have come to light only recently. For the 
values of A that arise in such compactifications, Chen et 
al. [H have shown that late-time 'eternal' acceleration is 
possible for open universes (see also Q). Andersson and 
Heinzle have further shown that 'eternal' acceleration is 
possible in the stronger sense of "at all times" (and not 
merely in the future of some initial onset of acceleration), 
and that there is a unique open universe cosmology with 
this property [l(|. Here we extend this result, in the form 
of a conjecture, to other values of A, including those that 
arise from flux compactifications. 
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II. PRELIMINARIES 

We consider a D-dimensional FLRW spacetime with 
metric 



s 2a^(r) / 2 (r)dT 2 + e 2/3 ¥ >(r) ds 2 



(III) 



where dE 2 , is the metric for a (D — l)-dimensional 
maximally symmetric space of constant curvature k = 
— 1,0, 1, and a, (3 are the D-dependent constants 



Q 



1 (D-l) 
2{D-2) ' 



0=l/y/2(D- !)(£>- 2). (II.2) 



Supposing that a and x are functions only of r, we find 
that the field equations reduce to equations of motion, 
and a constraint, for the variables ((p, ct, x) that are the 
Euler-Lagrange equations of the effective Lagrangian 



Leff = \ [r 1 (<P 



2 a 2 



-fe 2aip (k(3~ 



-20 V 



- 2Ae- ACT )] (II.3) 



where the overdot indicates differentiation with respect 
to t. We may take / to be any positive function. We 
may also make any choice of the positive constant A since 
its value is changed by a shift of a. It is convenient to 
choose 



f = e 



and to define 



2A = 1 



(II.4) 



(II.5) 



In this case, the equations of motion following from the 
effective lagrangian are equivalent to the coupled ordi- 
nary differential equations 



1 



1 



u = -auv + ^(A - fi)(u 2 + 1) + -fiv 2 + 



-av 2 H — Xuv + a 



Ilk 

2 r " ' 2/32 

\a-2f3ip 



2" 2af3 2 
together with the constraint 



X 2 = e-^ 



1 . l„A ff -2^ 
+ (3 2 



(11.6) 



(II.7) 



The x equation of motion is implied by differentiating 
the constraint. 

The spacetime metric can be put in the standard 
FLRW form 



ds 2 



-d^ + S 2 ^) d£ 



by defining 



S(t) = e^ (r) 



(II. 



(II.9) 



The condition for the universe to be expanding is dS/dt > 
0, which is equivalent to (p > 0. The condition for positive 
acceleration is d 2 S/dt 2 > 0, which is equivalent to 



v 2 < 2a 2 



(11.10) 



Using the constraint to eliminate the fc-dependent term, 
one finds that the acceleration condition is equivalent to 
the inequality 



< aj3 . 



(11.11) 



This has a simple physical interpretation because the left 
hand side is the kinetic energy of the scalar fields. 

We will now consider in turn the two special cases of 
(i) constant axion and (ii) k = 0. In these two cases the 
problem can be reduced to the analysis of a 2-dimensional 
autonomous dynamical system. For both cases, it is con- 
venient to define 



A c = 2V"/3 = 



\ h =2a = VD-l\ c 



D-2 1 

( IL12 ) 

These quantities have a general convention-independent 
significance, but take the above particular values in our 
conventions. 



III. CONSTANT AXION 



It is consistent to seek solutions of (|II.6[) and (|II.7|) by 

setting x = 0, in which case it is convenient to use the 
constraint to eliminate the fc-dependent terms from the 
equations of motion. The resulting equations of motion 
are 



-auv + 2^( u 



1 



-P(v 2 



and the constraint is 



2 2 
u — v 



2a 



-Xuv, (ILIA) 



1 _ ^_ p \<J-2Pv 

I3 2 



(m.2) 



Equations (|III.1|) define a two-dimensional autonomous 
system, analyzed by Halliwell for D = 4 [ll[. The gen- 
eralization to arbitrary D was given in In these ref- 
erences, the phase portraits were presented in the (u, v) 
plane. Here we summarize the results of a global phase 
plane analysis [l2j, paying particular attention to cos- 
mic acceleration; setting x — in pi. lip , we find that 
acceleration occurs when 



< 2a/3. 



(1TI.3) 



By using alternative phase plane coordinates one can 
study the nature of the trajectories at infinity, fixed 
points in particular. Actually, since trajectories at infin- 
ity are the same for positive and negative A, the phase- 
plane portrait for A > can be joined to the phase-plane 
portrait for A < at infinity. The resulting 'total' phase 
plane is topologically a sphere. The fixed points on the 
sphere have positions that are functions of A and they all 
are hyperbolic for finite (non-negative) A except when 

• (i) A = A c . The k = and k ^ fixed points 
in the (u, v) plane coincide at this 'critical' value 
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of A, which is a transcritical bifurcation point in 
the family of dynamical systems parametrized by A 
0. 

(ii) A = Aft. This was called the 'hypercritical' 
point in 4]. The k = fixed point goes to in- 
finity at this value of A. Globally, this value of A is 
another transcritical bifurcation [12j . 



where the second inequality is saturated for D = 9. 




FIG. 1: A c < A < A < Aft, and D < 9. The acceleration 
region is shaded. Regions A± correspond to open universes 
that are either expanding (A+) or contracting (A-). Region 
B corresponds to closed universes. 

Here we are concerned only with the hemisphere that 
contains the trajectories of the A > model, includ- 
ing the circle at infinity. This is illustrated for A > A c 
in Figs. 1&2, where the shaded region is the accelera- 
tion region. Fig. 1 applies for A < Aft and Fig. 2 for 
A > Xh- These two cases arise naturally in compacti- 
fications from general relativity in a higher dimension, 
with A < Aft arising from compactification on hyperbolic 
spaces and A > Aft arising from flux compactifications 
(see e.g. (HI). Trajectories in region A + correspond to 
expanding open (k = —1) universes, while those in region 
A_ correspond to contracting open universes. Region B 
trajectories correspond to closed (k = 1) universes, which 
may be contracting or expanding. Note that there are 
two fixed points that are simultaneously in A + and on 
the boundary of the acceleration region. The fixed point 
at infinity is a saddle and there is a separatrix trajectory 
that connects it to the other of these two fixed points, 
which is either a stable focus or a stable node depending 
on both A and D. This is a general feature for A > A c , 
which we assume in the following. As shown in this 
k = — 1 fixed point at finite (u,v) is a stable focus when 
both D < 9 and 



X> X = 



^(10 -£>)(£> -2) 
Otherwise it is a stable node. Note that 
A c < A < Aft 



(III.4) 



(III.5) 




FIG. 2: A > Aft, and D > 10. The acceleration region is 
shaded. Note the (separatrix) trajectory of region A+ that 
lies entirely within the acceleration region. 



A node is diffcomorphic to a focus and, consequently, 
they are often regarded as equivalent in texts on dynam- 
ical systems. However, the distinction is physically rel- 
evant in the present context because a diffeomorphism 
would alter the boundary of the acceleration region. In 
the approach to a focus on the boundary of this region a 
trajectory will ultimately cycle in and out of the region. 
In contrast, in the approach to a node on the boundary 
of the acceleration region, a generic trajectory will ul- 
timately be continuously inside or continuously outside 
the region. Moreover, for a node the possibility arises 
that the above-mentioned separatrix trajectory will lie 
entirely in the acceleration region; no other trajectory 
can have this property because the fixed point at infin- 
ity on the acceleration boundary is a saddle. A natural 
conjecture is that the separatrix trajectory lies entirely 
inside the acceleration region if and only if the 'internal' 
fixed point is a node rather than a focus. This is illus- 
trated for A < Aft in Fig. 1, and for A > Aft in Fig. 
2. 

This conjecture is consistent with recent results of An- 
dersson and Heinzle [io| . as we now explain. These au- 
thors considered open universes in D = 1 + m dimen- 
sions (m > 2) arising from compactification on an n- 
dimensional compact hyperbolic space. The equations 
that they consider are the same as those deduced from 
an effective £>-dimensional Einstein-dilaton theory with 
(see e.g. [lj|) 



2(m + n- 1) _ , 
A = a -, -t— = Xah ■ 

y nyra — 1) 



(n > 2). (III.6) 



Note that A c < Xah < Xh, so the fixed points are as 
shown in Fig. 1 rather than Fig. 2. The global phase 
space of Andersson and Heinzle is region A+ in Fig. 1. 
They proved that the separatrix trajectory lies entirely 
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within the acceleration region if and only if n+m+l > 10. 
This result is consistent with our conjecture (and proves 
it for some cases). To see this, we first note that the 
condition n + m+l > 10 is automatically satisfied if 
m > 9, but then D > 10 and it is also true that the 
'internal' fixed point is a node. For m < 8, we have a 
node only if A < A, but one can show that 



Xah < A +> 



2(9- 



m) 



n(9 — m) 



< 



(III.7) 



and hence that the 'internal' fixed point is a node if and 
only if n + m + 1 > 10. Thus, our conjecture is supported 
by the results of [l(| but predicts the existence of an 
eternally accelerating trajectory for D < 8 as long as 
A > A, and also for A > A/ t when D > 10, as shown in 
Fig. 2. As mentioned above, this latter case is applicable 
to flux compactifications. 

It should be noted that eternal acceleration, in the 
above sense, contradicts observations, which also favour 
a flat universe, so eternally accelerating k = — 1 universes 
are of little astronomical interest. For the remainder of 
the paper we consider only flat universes, and we seek 
solutions that undergo multiple periods of acceleration. 



IV. VARYING AXION: FLAT COSMOLOGIES 

Another case in which the dynamics is reducible to 
that of a two-dimensional autonomous dynamical system 
is that of spatially flat cosmologies, but now allowing 
for arbitrary motion on the dilaton-axion target space. 
Setting k = in (|II.6p we have the equations 

u = —auv + — (A — n) (u 2 + 1) + 2^ v2 ' 



-av 2 + —Auv + a , 



(IV.l) 



which define a 2-dimensional autonomous dynamical sys- 
tem. Setting k — in (|II.7[) we have the constraint 

X 2 = e-» a [v 2 - u 2 - 1] , (IV.2) 

which restricts the physical phase space by the inequality 

v 2 - u 2 - 1 > . (IV.3) 

Note that this is a non-compact region of the (u, v) plane. 
Before proceeding, it is convenient to introduce new vari- 
ables that will allow a compactification of the phase space 
and thus facilitate the analysis of any fixed points at in- 
finity. We choose these 'global' variables to be 1 



y = v 



u = xy 



v = y 



The constraint now implies that x 2 + y 2 < 1, so that the 
allowed region of the (u, v) phase plane has been mapped 
to the unit disk in the (x, y) plane. If we define a new 
independent variable f such that 



\v\dr = dr . 



(IV.5) 



then the equations that define the dynamical system in 
the (x, y) variables are 



dx 
~dJ 
dy 
dr 



sgnl//) { y 2 
1 



(A — fi) — ax 



\y\< a - -Xx - ay' 



(IV.6) 



The sign of y is fixed for any given trajectory, so these 
equations define two independent dynamical systems, one 
for y > and another for y < 0. As the condition for 
expansion in the new variables is y > 0, trajectories with 
y > correspond to expanding universes while those with 
y < correspond to contracting universes, and time re- 
versal takes one set of trajectories into the other 2 . The 
condition for positive acceleration (which is unchanged 
by time reversal) is (jll.lOj) . but the equation of motion 
for v can be used to simplify this to v 2 < 2a 2 , which is 
equivalent to 



y 



2 > 2/X 2 



(IV.7) 



Fixed points associated with infinite values of u and v 
have the following (x, y) coordinates, and eigenvalues: 



( tl.lh : {a+^+fi 



(IV. 



The nature of the fixed points at infinity therefore de- 
pends on the values of A and /i. Excluding /i = and 
A ^ Ad, which yield non- hyperbolic fixed points, we have: 

• A < Ah, li > 0: (1,0) is a saddle, and (—1,0) is an 
unstable node. 

• A < Ah, fi < 0: (1,0) is an unstable node, and 
(—1,0) is a saddle. 

• A > A/i, /i > 0: (1,0) is a stable node, and (—1,0) 
is a saddle. 

• A > Ah, /J < 0: both (1,0) and (—1,0) are saddles. 

For A < Xh there is a pair of fixed points, not at infinity, 
on the boundary of the allowed region, and hence on the 
unit circle in the (a;, ?/)-plane; one has y > and the 
other y < 0. The (x, y) coordinates, and eigenvalues, of 
the y > fixed point, are 



(IV.4) 




A 2 -A 2 A( Mc -/■!) • 
2Xh Xh 



, (IV.9) 



The phase portrait is invariant (up to direction of the trajecto- 
ries) under reflection about the x-axis so the variables (x, y) with 
y = y 2 constitute another possible choice of global variables; this 
choice was made in but our choice makes the issues of expan- 
sion and acceleration clearer. 



2 Each of these two dynamical systems can be compactified by the 
addition of a y = trajectory but the two y = trajectories 
cannot be identified because their directions are opposed. The 
full compactified phase space is therefore two disjoint hcmicircles, 
as shown in Figures 3 and 4. 
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respectively, where 



A. A > A;,: Recurrent Acceleration 



= j (A 2 - A 2 ) 



(IV.10) 



This type of fixed point is a stable node for /i > [i c 
and a saddle for [i < fi c . It is non- hyperbolic on the 
curve /i = /i c (A), but we will leave consideration of these 
cases to another article [l2j]. We note, however, that the 
Freedman-Schwarz supergravity model mentioned in the 
introduction has \i = fi c , with A = A c , so this is an ex- 
ample of a model with non- hyperbolic fixed points. 

Finally, if — A) > then there is a fixed point for 
y > in the interior of the unit circle, and another for 
y < 0. The (x, y) coordinates, and eigenvalues, of the 
y > fixed point are 

X h 



JJ_ 

x — fi ' y fi — x 

respectively, where 



{E+,E_} 



(IV. 11) 



2(/i - A) 



(«±Va7^) (IV. 12) 



with 



A(A, m) = a 2 fi c + 2A(/i - ^ c ) 2 - 7a 2 (n - fi c ) . (IV.13) 

Given that this 'interior' fixed point exists, which it will 
do only if — A) > 0, then it lies within the allowed 
region only if fi < fi c . Together, these two conditions arc 
equivalent to 

fj, < inf[0, a« c ] . (IV. 14) 

As this implies that [i < 0, the fixed point will be a stable 
focus whenever A > 0. 

At the 'interior' fixed point we may integrate the equa- 
tions & = u and ip — v to get the solution 



2a 



V7*G"-A) 



( T ~ T o) , 



X 



± 



y/\(Hc - /i) 

a/j, 



(t-tto) 



'/' - 



A 



Xo • 



(r - T ) + ifio . 



(IV.15) 



By an appropriate choice of ^o, the metric can be brought 
to the standard FLRW form ([11.81 ) with scale factor 

A 



S = t s , 
where the FLRW time is 
t oc exp 



5 = 



A* 



X(D-1) 



Xa(r - T ) 



(IV.16) 



(IV. 17) 



.\W - A) 

This fixed point lies inside the acceleration region if 5 > 
1, which is equivalent to 

n/X<-2/X 2 c . (IV.18) 

Note that S — > oo as // — > — oo, so that the acceleration 
at the fixed point can be arbitrarily large. 

The implications of this 'interior' fixed point solution 
for the other k = trajectories is what will mostly con- 
cern us for the remainder of this article. The phase por- 
trait depends on whether there is also a fixed point on 
the unit circle, which there will be only if A < A/j. We 
will consider in turn the two cases A > A^ and A < A/,. 



For A > Xh, we have \i c > 0, so that (|IV.14[) becomes 
/i < 0. Since /i < fj, c , equation (jIV.13j) implies that 
A > 0. The 'interior' fixed point is therefore a stable 
focus. For /i < the fixed points at infinity are both 
saddles, so the 'trajectory' at infinity, y = 0, and the 
'boundary' trajectory x 2 + y 2 = 1, with y > 0, constitute 
a heteroclinic cycle enclosing the stable focus. All other 
k = trajectories with y > must approach the stable 
focus at late times. Going backwards in time, each such 
trajectory describes a quasi-cyclic early universe that ap- 
proaches the heteroclinic cycle more closely in each cycle, 
as shown in Figs. 3&4, where (in each of these figures) 
the lower-hemisphere corresponds to the time-reversed 
trajectories with y < 0. This behaviour was noted in 
2] for the D — 4 model with A = — /i = 2, except that 
they chose to consider the case for which the heteroclinic 
cycle is approached at late times; this happens for y < 
but then the Einstein-frame universe is contracting. As 
we insist on expansion in Einstein frame, we must choose 
y > 0, and then the approach to the heteroclinic cycle 
is at early times, going backwards in time. Of particular 
interest here are the implications for cosmic acceleration. 
The "acceleration region" is the shaded region in Figs. 
3&4. 




FIG. 3: A > Xh, /J- < /J-c, n/X > —2/Xc. Trajectories within 
the upper (lower) hemicircle represent expanding (contract- 
ing) universes. In each hemicircle, there is a focus outside the 
(shaded) acceleration region. Note the recurrent acceleration 
of generic expanding universes. 



The focus (which is stable for y > and unstable for 
y < 0) lies inside the acceleration region if the ratio fi/X 
satisfies pV.18|) . Since the quasi-cyclic trajectories only 
arise in the range of parameters for which fi c > 0, any 
value of n that satisfies ([IV.18[) will also be such that 
fi < [i c . Thus, there exists a range of parameters (fj,, A) 
such that there is a stable focus inside the region of accel- 
erated expansion, and another range that places it inside 
the region of decelerated expansion. These two possibil- 



ities are illustrated in Figs. 3&4 respectively, and the 
respective parameter ranges are shown in Fig. 5. 




FIG. 4: A > Xh, fJ- < He, /U/A > — 2/Ac- For each hemicircle, 
there is now a focus inside the acceleration region, so that a 
generic expanding universe undergoes a continuous late-time 
acceleration, preceeded by early-time recurrent acceleration 
and deceleration. 



For generic initial conditions set for very small scale 
factor (at which either the approximations of general 
relativity break down or the k = assumption fails), 
the universe expands with a series of rapid oscillations 
between acceleration and deceleration. Potentially, this 
might have the same effect as an early-time inflationary 
epoch if the number of cycles of acceleration is sufficiently 
large, although it must be remembered that the universe 
under discussion is always flat, by hypothesis, so this 
"quasi-inflationary" epoch could not be used to explain 
why our Universe is flat. At late times the oscillation 
between acceleration and deceleration gives way to con- 
tinuous deceleration, as illustrated in Fig. 3, or continuous 
acceleration, as illustrated in Fig. 4. For the special case 
that {j, = — 2A/A^, the stable focus lies on the bound- 
ary of the acceleration region and the oscillation between 
acceleration and deceleration continues forever. 

Finally, we recall that these results depend strongly 
on the fact that we have chosen the Einstein conformal 
frame. If we define a conformally rescaled metric 
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(IV.19) 



for arbitrary constant r, then we get a Lagrangian density 
of the form 



C = e r V-dett/E + 
and the FLRW scale factor becomes 



(IV.20) 



(IV.21) 




FIG. 5: (A, /x)-plane. If the pair of parameters (fi, A) is in- 
side the semi-infinite grey region, there exists a fixed point 
inside the unit circle. If (^, A) is in region I (light grey), this 
fixed point is outside the acceleration region. If (fi, A) falls 
into region 77 (dark grey) it is located inside the acceleration 
region. 



If rXh < 1 then all y > universes are still eternally 
expanding but for rXh > 1 there is a late-time oscillation 
between expansion and contraction, as observed in [2j for 
D = 4 and r = 1. 



B. A < X h 

For A < Xh the phase portrait is qualitatively different 
because of the 'finite' fixed point on the boundary of the 
allowed region. This fixed point (which is the same as the 
k = fixed point occurring for A < Xh in the analysis of 
the x = case of section 3) lies within the acceleration 
region if A < A c and outside it if A > A c . We shall 
consider only the A > A c case. 

For A < Xh we have fx c < 0, and hence, from (|IV.14jl . 
there exists an 'interior' fixed point when /j, < fi c . This 
fixed point will be a focus if 



1 

16A 



X h y / 81A2 _ 32A 2 + %X\ - 16A 2 



The universe is now expanding in the new conformal 
frame if v > rXhU or, equivalently for y > 0, rXhX < 1. 



(IV.22) 

Note that /i_ < /i c for A < Xh, so that there is a range 
of values of fi for which this condition is satisfied. This 
focus is again stable for y > and it again lies within 
the acceleration region if /i < ~2X/X 2 C . This condition is 
satisfied for A sufficiently small, in particular for A > A c 
provided that A— A c is sufficiently small. For larger values 
of A < A/j, the focus must lie in the deceleration region. 
Thus, given that the late-time behaviour of generic tra- 
jectories is governed by an 'interior' focus, this focus may 
lie inside or outside the acceleration region, depending on 
the values of the parameters (A, /i). The two possibilities 
are illustrated in Figs 6 &7. Inspection of these phase 
portraits shows that almost all trajectories begin at the 
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FIG. 6: A > \h, A* < Mc, A*/A < — 2/Aj?. Trajectories within 
the upper (lower) hemicircle represent expanding (contract- 
ing) universes. Within each hemicircle there is a focus outside 
the acceleration region. 



unstable node at (x,y) = (1,0) and end at the 'interior' 
focus (the separatrix trajectory connecting the two fixed 
points is the unique exception for which \ is n °t always 
zero). Thus, given the existence of the 'interior' focus, 
the late time behaviour for A c < A < A/j is similar to the 
late-time behaviour for A > Aft. 




FIG. 7: Ac < A < X h , n < fi c , fj,/X < -2/\ 2 c . Within each 
hemicircle there is now a focus inside the acceleration region. 



V. DISCUSSION 

We have presented a study of the homogeneous 
and isotropic cosmological models that arise in a two- 
parameter class of models of gravity coupled to dila- 
ton and axion fields in arbitrary spacetime dimension D; 



the two parameters (A, fi) determine the dilaton coupling 
constant (and hence the 'strength' of an 'exponential' po- 
tential) and the radius of curvature of the 2-dimensional 
hyperbolic target space (which we assumed to be finite). 
It is well-known that all possible cosmological solutions 
of the zero-axion model can be understood in terms of 
trajectories of a two-dimensional autonomous dynami- 
cal system, but we have explored some new features that 
arise from consideration of the global nature of the phase 
space of this system. It is less well-known that a similar 
reduction to a two-dimensional autonomous dynamical 
system is possible for the generic dilaton-axion model 
provided that the restriction is made to flat cosmologies; 
a particular D = 4 dilaton-axion model was analysed by 
this method in and we have extended this analysis to 
the generic dilaton-axion model. In addition, we have ex- 
plored here the consequences for Einstein frame cosmol- 
ogy, which are quite different from those in the 'string- 
inspired' frame of 2]. In particular, all flat Einstein- 
frame cosmologies are either eternally expanding or eter- 
nally contracting, but may undergo one or more periods 
of cosmic acceleration. Our main result is that there is 
a large range of the parameters (A, fi) for which almost 
all flat expanding universes undergo an epoch of oscil- 
lation between acceleration and deceleration, with cither 
continuous late-time acceleration or continuous late-time 
deceleration depending on the value of the ratio X//i (and 
eternal oscillation at a critical value of this ratio). 

The early-time behaviour depends strongly on whether 
A is greater or less than the 'hypercritical' value A^. For 
A > Xh, the global phase space is topologically equiva- 
lent to the one given in Q, in which almost all trajec- 
tories are asymptotic in one direction to a heteroclinic 
cycle and in the other direction to a focus within this 
cycle. The focus can be stable or unstable according to 
the chosen direction of time. In a non-Einstein frame, 
such as that considered in 0], for which the universe 
undergoes quasi-cyclic periods of expansion and contrac- 
tion, one may suppose that the focus is unstable and 
the the heteroclinic cycle is approached at late times. In 
the Einstein frame, for which the universe is either eter- 
nally expanding or eternally contracting, the heteroclinic 
cycle is approached at late times only by a contracting 
universe. If we assume expansion, then the heteroclinic 
cycle dominates the early-time behaviour, which exhibits 
oscillation between acceleration and deceleration. This 
oscillation continues until the focus is approached suffi- 
ciently closely that its behaviour, of continuous acceler- 
ation or continuous deceleration, takes over. This is the 
origin of the epoch of oscillation between cosmic acceler- 
ation and deceleration when A > Xh ■ For A < A^ there is 
no heteroclinic cycle and the early-time behaviour of ex- 
panding universes is quite different, but the medium-time 
behaviour may again exhibit oscillation between cosmic 
acceleration and deceleration for values of the parame- 
ters that put the focus sufficiently close to the borderline 
between acceleration and deceleration. 

The fact that there are models for which almost all flat 
universes undergo cycles of cosmic acceleration and de- 
celeration is of little use for inflationary purposes because 
one aim of inflation is to explain why the Universe is flat, 
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and for that purpose one needs to consider universes that 
are not flat. In this case, however, the cosmological evo- 
lution cannot be reduced to the study of a 2-dimensional 
phase-plane. It would be of interest to use other methods 
to learn more about generic non-flat cosmologies (with 
non-constant axion) but at present we must assume that 
there exist some non-flat cosmologies that become ap- 
proximately flat at late-times, in which case their late- 
time behaviour will determine the 'initial' conditions for 
a flat cosmology. Thus, the early-time behaviour of the 
flat cosmologies studied in this paper is probably not 
physically relevant. The late-time behaviour may not 
be relevant either, depending on whether an almost flat 
universe becomes more or less flat as time progresses. It 
is plausible that flatness is a stable feature of universes 
that exhibit continuous late-time acceleration, and an 
unstable feature of those that exhibit continuous late- 



time deceleration. However, the medium-time recurrent 
acceleration that we have found to be typical of many 
dilaton-axion models could be relevant to the observed 
cosmic acceleratiion of our Universe. At the very least, 
it shows that current observations of cosmic acceleration 
provide no evidence that this acceleration will continue 
forever. With many scalar fields, the future is uncertain! 
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